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We study by the perturbative Functional Renormalization Group (FRG) the Random Field and 
Random Anisotropy 0{N) models near d = 4, the lower critical dimension of ferromagnetism. The 
long-distance physics is controlled by zero-temperature fixed points at which the renormalized ef- 
fective action is nonanalytic. We obtain the beta functions at 2-loop order, showing that despite 
the nonanalytic character of the renormahzed effective action, the theory is perturbatively renor- 
malizable at this order. The physical results obtained at 2-loop level, most notably concerning the 
breakdown of dimensional reduction at the critical point and the stability of quasi-long range order 
in d < 4, are shown to fit into the picture predicted by our recent nonperturbative FRG approach. 

PACS numbers: ll.lO.Hi, 75.40.Cx 



I. INTRODUCTION 

Despite decades of intensive investigation the effect of 
weak quenched disorder on the long-distance physics of 
many-body systems remains in part an unsettled prob- 
lem. This is the case for the class of models in which 
iV-component classical variables with 0{N) symmetric 
interactions are coupled to a random field. Depending 
on whether the coupling is linear or bilinear, the mod- 
els belong to the 'random field' (RF) or the 'random 
anisotropy' (RA) subclasses. Such models with = 1, 2, 
or 3 are relevant to describe a variety of systems encoun- 
tered in condensed matter physics or physical chemistry. 
To name a few, one can mention dilute antiferromagnets 
in a uniform magnetic fieldji critical fluids and binary 
mixtures in aerogels (both systems being modelled by the 
A^ = 1 RF Ising model) (2*2*4 vortex phases in disordered 
type-II superconductors (described in terms of an elas- 
tic glass model whose simplest version is the N ^ 2 RF 
XY model) f^iSiSiS amorphous magnets, such as alloys of 
rare-earth compounds^Siii and nematic liquid crystals in 
disordered porous media (described by A^ = 2 or A^ = 3 
RA models) 

On the theoretical side the main questions raised about 
the equilibrium behavior of such systems concern the 
nature and the characteristics of the phases and of the 
phase transitions. It has been shown by both heuristic 
and rigorous methodsi^ii^iiSiiSiii that the lower critical 
dimension below which no long-range order is possible 
is die = 2 for the RFIM and die = 4 for RF models 
with a continuous symmetry {0{N) with A^ > 1). The 
same conclusion applies to RA models with the restric- 
tion that only ferromagnetic (to use a magnetic termi- 
nology) long-range order is forbidden below c?ie = 4; an- 
other type of long-range order associated to a spin-glass 
phase is still possible>ii Here, we only consider RA mod- 
els with isotropic distributions of the random anisotropics 
and with A^ > 1; for anisotropic distributions, long-range 
ferromagnetic ordering may still occur below die = 4, 
whereas RA makes no real sense for A^ = 1, the model 



reducing then either to the random temperature Ising 
model or to the pure Ising model depending on the de- 
tails of the effective hamiltonianii 

Two central issues remain under active debate. The 
first one is about the so-called 'dimensional reduction' 
property. Standard perturbation theory predicts to all 
orders that the critical behavior of an 0{N) model in 
the presence of RF is the same as that of the pure 
model, i.e. with no RF, in two dimensions less^ 
The same applies to the RAO (TV) M with A^ > 1 
near the paramagnetic-ferromagnetic transitionii^ Di- 
mensional reduction is known however to break down, 
its most striking failure being the prediction of a lower 
critical dimension die = 3 for the RFIM in contradiction 
with the exact result (see above). A proper description 
of the long-distance behavior of RF and RA models must 
thus provide a way out of the dimensional reduction. 

The second issue concerns the phase diagram of the RF 
and RA models with a continuous symmetry (A^ > 1) in 
dimensions below d — 4, which of course are relevant 
to the physical situations. If long-range ferromagnetism 
is forbidden, quasi- long-range order (QLRO), namely a 
phase characterized by no magnetization and a power-law 
decrease of the correlation functions at large distances, 
may still exist i^*SiLSi2i2a It has been shown that QLRO 
is absent for A^ > 3 in the presence of RF and for A^ > 
10 in the presence of RAjSS yet it has been argued that 
QLRO is present for A^ = 2 in d = 3, in which case 
it corresponds to the 'Bragg glass' phase predicted for 
vortices in disordered type-II superconductorsi^iSiS 

We have recently proposed a coherent resolution of 
those issues based on a non-perturbative (NP) functional 
renormalization group (FRG) treatment iSiiS2 This ap- 
proach has allowed us to provide a unified picture of fer- 
romagnetism, QLRO and criticality in RF models in the 
whole (A^, d) diagram as well as a way to escape dimen- 
sional reduction. 

The main finding a^^i^^ can be summarized on the phase 
diagram of the RFO(Af)M shown in Figure Q In re- 
gion III, there are no phase transitions and the system 
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FIG. 1: Nonperturbative FRG prediction for the (A'^, d) phase 
diagram of the RFO(Af)M. See text for comments. 



is always disordered (paramagnetic). In regions I and 
IV, there is a second-order paramagnetic to ferromag- 
netic transition and in region II, a second-order transition 
between paramagnetic and QLRO phases. In all cases, 
the critical behavior is controlled by a zero-temperature 
fixed point at which temperature is formally irrelevant. 
At this fixed point the renormalized effective action is a 
nonanalytic function of its arguments (the order parame- 
ter fields). Although present, the nonanalyticity is weak 
enough in region IV to let the critical exponents take their 
dimensional reduction value (corrections to scaling may 
nonetheless differ from the dimensional reduction predic- 
tions). In regions I and II the nonanalyticity takes the 
form of a cusp in the renormalized second cumulant of 
the random field, which leads to a complete breakdown of 
dimensional reduction. Finally, the whole QLRO phase 
in region II is also controlled by a zero-temperature fixed 
point characterized by a cusp. 

There is undoubtedly room for improving the quan- 
titative predictions of our NP-FRG theory, in terms of 
both the number of observables studied and, more im- 
portantly, of the accuracy of the (necessary) approxima- 
tions involved. (In addition, the NP-FRG study of the 
RAO(A^)M has not yet been completed.) The robust- 
ness of the proposed scenario may however be tested by 
considering a perturbative FRG treatment of the models 
near d = 4. Such a perturbative FRG has been pio- 
neered by D. FisheaiSiSi and widely used to study the 
statics and the depinning of elastic systems pinned by 
quenched disorderi ^i^i^i^^^i^^ 

At one-loop level, the flow equation for the renormal- 
ized second cumulant of the disorder, first derived by D. 
Fisher— for the RF and RA 0{N) models, has been stud- 
ied by several authors. The results fit into the 
diagram displayed in Figure Q which should come as no 
surprise: the flow equations obtained in our NP-FRG ap- 
proach exactly reproduce the 1-loop result near d = Am^ 
Below, we give a survey of the behavior of the RFO(A^)M 
at 1-loop level in d = 4 -I- e, including some new results, 
as well as a study of the related RAO(A^)M. 



To go beyond this first step, one must consider the 
next order in the loop expansion. However, the tech- 
nical difficulties are now much more involved than at 
the 1-loop level. On top of the rapidly increasing num- 
ber of diagrams, diagrams which in the present case are 
functionals, the nonanalytic character of the renormal- 
ized effective action at T = leads to the appearance of 
'anomalous' terms in the diagramatics, whose evaluation 
is a priori ambiguous. A resolution of the problem has 
been proposed for the simpler case of disordered elastic 
systems by Le Doussal, Wiese and coworkers— 

A preliminary account of this work has been published 
in R ef. 1211 An independent calculation has appeared in 
Ref.il 

The rest of the paper is organized as follows. In sec- 
tion m we present the RF and RA 0{N) models and 
their nonlinear sigma versions appropriate to describe the 
long-distance physics near the lower critical dimension of 
ferromagnetism, d — 4. We outline the perturbative FRG 
framework and the way to extract scaling behavior and 
critical exponents. Section lllll is devoted to an analy- 
sis of the 1-loop FRG equations at T = in d = 4 -I- e, 
discussing the fixed points and their stability and con- 
trasting the RF and RA cases. In section Hvl we derive 
the FRG beta functions at 2-loop order in T = 0. We 
present the diagrammatic representation and the method 
used to handle the apparent ambiguities appearing in 
the formulation due to the nonanalytic character of the 
renormalized dimensionless effective action. Proceeding 
in this way we obtain a well-defined renormalized the- 
ory at 2-loop order. The physical results obtained from 
solving the 2-loop FRG equations are discussed in sec- 
tion and we stress the new features appearing at this 
order. Several technical aspects of the work are deferred 
to Appendices. 

II. MODELS AND FRAMEWORK 
A. Models 

We consider the 0{N) model in the presence of RF 
or RA near d — A. We stress again that d — 4 is the 
lower critical dimension for A'^ > 1 and for the param- 
agnetic to ferromagnetic transition. In a manner similar 
to that developed for the pure model at low temperature 
near d = 2, the long-distance physics for weak disor- 
der (which, we recall, takes here the role played by low 
temperature in the pure model, temperature being now 
irrelevant and eventually set to zero) can be described 
in a field-theoretical setting by a nonlinear sigma model 
with effective Hamiltonian 

n[S] = J d'^x l(VS(x))2-^(/,/(x)+i/u')5^(x) 
-^r'^(x)5*(x)5^"(x) 

ij 

(1) 
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where the iV-component spins S satisfy the fixed-length 
constraint, S(x)^ = 1, and H = Hu is a uniform exter- 
nal magnetic field; h(x) is a random magnetic field and 
t(x) a second-rank random anisotropy tensor, both with 
gaussian distributions characterized by zero means and 
variances given by: 



where the 1-replica and 2-replica parts of the action read 



/i'(x)/iJ(y) = A5y 5(x-y) 



T'J (x)Tfc'(y) = -f (S.kSji + 6aS,k) ^(x - y). 



(2) 



(3) 



Higher-order random anisotropics could be included as 
well. They will indeed be generated in the perturbation 
expansion and the renormalization group flow.^^ How- 
ever, starting with only a second-rank (or more gener- 
ally an even-rank) random anisotropy, only even-rank 
anisotropics will be generated. In what follows we will 
therefore use the acronym RA to characterize models 
with even-rank random anisotropics. 
From the associated partition function. 



VS S{S^ - 1) exp 



(4) 



one can obtain the free energy by averaging the loga- 
rithm of Z over the quenched disorder. This is more 
conveniently performed by introducing replicas Sq(x), 
a = l,...,n, which leads after explicitly performing the 
average over the disorder to the following 'replicated' ef- 
fective Hamiltonian 



H„[{SJ]= f rf'^x^i(VS(x))^-5]i/u.S(x) 

a a 
-^E^O (Sa(x)-Sfc(x)) 



(5) 



ab 



with 



(6) 



and —l<z<+l. The fluctuations around a fully or- 
dered state in which spins in all replicas align in the same 
direction are as usual handled by splitting the replica field 
Sa(x) into a component collinear to the external field 
(and to the magnetization), Sa(x) = Sa(x).u, and A^ — 1 
components orthogonal to it, na(x) = Sa(x) — Ea(x)u. 
By using the relation between I]a(x) and na(x) imposed 
by the unit- length constraint, the replicated partition 
function can be expressed as a functional integral over 
the {N — l)-component replica fields Ila{^), 

•' a 

expf-^5i[n„] + i^52[na,n,,] + --- J , 

\ a ab / 

(7) 



5i[nj 



2(1 -n2) 



H 



(8) 



(9) 



and the dots denote terms, such as those produced by 
the Jacobian of the transformation from the S^'s to the 
Ila's and possible contributions involving more than two 
replicas, that either do not contribute to the perturbation 
expansion in the T = limit or turn out to be irrelevant 
within conventional power counting. ^^'^^ 

From the logarithm of the partition function, Eq. (Tj), 
one can obtain, by a Legendre transform with respect 
to external sources coupled to the {N — l)-component 
replica fields na(x), the effective action r„ [{IIq}] which 
is the generating functional of the onc-particle irreducible 
vertices for the Ila fields and from which all equilibrium 
observables can be derived. (The subscript n will be 
dropped in the following.) 



B. Perturbation theory and renormalization 

We proceed by calculating the effective action F [{Ha}] 
perturbatively in powers of the disorder correlator Rq, 
keeping only terms that do not vanish in the limit T = 0. 
The results so obtained would however be singular, show- 
ing the standard ultra-violet divergences as e = d— 4 goes 
to zero. For instance, if we use the dimensional regular- 
ization as a regularization scheme, the 1-loop calculation 
brings in terms proportional to 1/e. To cure this problem, 
it is necessary to renormalize the theory by introducing in 
the effective Hamiltonian 'counterterms' that are chosen 
to keep the physical quantities finite. 

Expressed in terms of dimensionless renormalized 
quantities at an arbitrary momentum scale fi, the 1- and 
2-replica parts of the action read: 



Si [iTa] 



d-2 y 



2ZTt 



(tTq • VTTa)) 



L 1 - Znni 

- 2Z^'^''ZTh^l-Znnl} 



(10) 



ZZlrpt J \ ^^^-^ 



^JT^~Z^^J 1 - ZnTT, 
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where the dimensionless renormahzed quantities are de- 
fined as 



(12a) 
(12b) 

(12c) 
(12d) 



the renormahzed and the bare theories is simply 



n = V^nTT 

1 

H 



and Zji[z) is a functional of the renormahzed dimension- 
less disorder correlator R{z), with its leading term equal 
to R{z). The two renormalization constants Zt and Zu 
and the renormalization function Zfj{z) are chosen so 
that the loopwise perturbative expansion of the effective 
action remains finite. (We work in the minimal subtrac- 
tion scheme in which the counterterms contain only the 
singular parts necessary to make the physical quantities 
finite.) In practice we compute the 2-point proper vertex 
associated with the 1-replica part of the effective action, 

(2) 

T\ '{q), and the 2-rephca part of the effective action, 
both being evaluated for uniform configurations of the ttq 
fields. 

The perturbative expansion is organized about the free 
theory formed by the quadratic part of the 1-replica ac- 
tion. The associated free propagators are expressed in 
terms of the bare quantities as follows: 



g2 + H/Y. 



(13) 



where S = 1 — 11^. Note that following Brezin and 
Zinn-Justin^^ we keep an external magnetic field H which 
allows to regularize the infrared divergences by giving a 
mass to the Goldstone modes. Aside from this term, 
the action in Eqs. (10,11) is 0{N) invariant. The loop 
expansion can be graphically expressed in terms of 1- 
particle irreducible Feynman diagrams with vertices com- 
ing from both the non-quadratic piece of the 1-replica 
action and from the 2-replica action. 

A difficulty of the present problem lies in the functional 
character of the expansion, the 2-replica vertices involv- 
ing the whole function R{z) and its derivatives. This is 
somewhat similar to the treatment of disordered elastic 
systemSfSi, with however the additional complication that 
the 1-replica part is now nontrivial and gets renormahzed 
in a manner that couples to the renormalization of the 
disorder. The details of the calculation as well as the 
method to handle possibly anomalous terms appearing 
at 2-loop level when the renormahzed correlator of the 
disorder is nonanalytic will be presented in Section Hvl 



C. FRG equations, critical exponents and 
correlation functions 

For the 1-replica, 2-point proper vertex and for the 
2-replica effective action (when both evaluated for uni- 
form configurations of the fields), the relation between 



t[%; tt, t, h, R) = Zu rYi,{q; H, T, H, i?o) (14) 



.(2) 



T2^f_,{z, t, h, R) = r2,B(2, r, H, Ro) 



(15) 



where B denotes the bare theory. 

The RG flow equations then result from the invariance 
of the bare theory under a change of the momentum scale 
/i, when T, H and Rq are held fixed. Actually, we are only 
interested in the situation of zero temperature (T = 0) 
and zero external field {H — 0). We introduce 



and 



Cn = ^J,^^log ZuIro 
Ct = fJ.d log Zt\ro 



Pb{z) = -tJ,d^R{z)\R„, 



(16) 
(17) 



(18) 



where we have implicitly set T = H = 0. As an illus- 
tration, the flow of the 1-replica proper vertex, tr\ ^{q), 
when H = h = 0, T = t = 0, U — TT = is derived as 



+ (2 - d + Ct - Cn)- 



(irS(9)) = o 



(19) 



where the long-distance physics is now obtained when 
^ ^ 0. 

The scaling behavior and the critical exponents of the 
physical quantities can be obtained from the fixed-point 
solutions and the properties of the flow near the fixed 
points. In particular, the exponents 77 and fj that char- 
acterize the power-law decay of the 2-point correlation 
functions at the critical point for small q, 



(S(-q) ■ S(q)) - (S(-q)) • (S(q)) ~ q 



(S(-q)) . (S(q)) - (S(-q)) • (S(q)) 



-(i-v) 



(20) 



(21) 



and the exponent 9 associated with the temperature, t ~ 
/i^T, are given by 



V ~ Cn* — Ct* 

^ = 4 - + Cn* 

e = d-2-CT,=2-rj 



(22) 
(23) 
(24) 



where Cn* and Ct* are the fixed-point values of Eqs. 
(|16I17|I . Provided 6* > 0, the fixed point indeed occurs at 
zero renormahzed temperature. 

Before closing this section, it is worth recalling an 
inequality for the correlation functions in the present 
models, which turns into an inequality between criti- 
cal exponents. In the RF case, the result is due to 
Schwartz and Sofferf2£ who have proven that the q 
Fourier component of the 'connected' pair correlation 
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function, (S(— q).S(q)) — (S(— q)) . (S(q)), is always less 
than the square root of the q com ponent of the 'd is- 
connected' pair correlation function, (S(— q)) . (S(q)) — 
(S(— q)).(S(q)), up to an irrelevant multiphcative con- 
stant. As a consequence, one must have fj < 2r]. 

The RA case is different and has been considered by 
Feldman.^*^ In this model indeed, the randomness couples 
to a composite field that is bilinear in the spin variables 
(see Eq. As a consequence, the inequality now ap- 

plies to the connected and disconnected correlation func- 
tions of the composite (bilinear) field. One can define 
new critical exponents, r/2 and 7/2, for those correlation 
functions, 

(m(-q) . m(q)) (m(-q)) • (m(q)) ~ q-^^-^^'' (25) 



(m(-q)) . (m(q)) - (m(-q)) • (m(q)) ^ g-^^-,,) (36) 

where m*(x) = S'*(x)^ — (1/A^). The inequality between 
the correlation functions then imposes that 772 < 2r/2. 
However, the exponents rj and fj are no longer constrained 
by the usual Schwartz-Soffer inequality. The exponents 
772 and 772 are also expressable in terms of fixed-point 
quantities. 

In the following we first analyze the FRG equations 
obtained at 1-loop order near d — A. 



]^ ..19.20 

V = CKil) (28) 

fj^-e+{N-l)CK{l) (29) 

772 = (iV-f 2)Ci?t(l) (30) 

fj2 = -e + 2NCKil) (31) 

where the star indicates a fixed-point solution. 

For studying the fixed points and their stability it is 
convenient to introduce R{z) ~ {C/e)R{z) (the renor- 
malized disorder is of order e at the putative fixed points) 
and to consider the beta function for its derivative, 

e-i/3j^,(z) = -R'iz) + zR'{zf + {N - 3)R'{z)R'{l) 
+ {N~3 + 4z2) R'{z)R"{z) -{N + l)zR'{l)R"{z) 

- z (1 - z^) R'{z)R"'{z) + (1 - z^) R\l)R"'{z) 

- 3z (1 - z") R"{zf + (1 - z^Yr"{z)R"'{z) 

(32) 

It is illustrative to write down the 1-loop beta functions 
for the first derivatives R'{z — 1) and R"{z — 1), assum- 
ing that R{z) is at least twice continuously differentiable 
around z = 1. {z = \ corresponds to the situation where 
the spins in the two considered replicas become equal, 
Sq Sfe.) The expressions are 

e-'l3~^^^~R'{l) + {N^2)R\lf (33) 



III. ANALYSIS OF THE 1-LOOP FRG 
EQUATIONS IN d = 4 + e 

A. 1-loop beta function in d = 4 4- e 

The beta function for the renormalized correlator of 
the disorder R{z) at zero temperature has been obtained 
by Fisher at the 1-loop level in c? = 4 + e.iS. it reads 

I3r{z) = ~fidf,R{z) ^ -eR{z) + c(2(A^ - 2)R{z)R'{l) 
+ i(7V - 2 + z^)R'{zf - z(l - z^)R'{z)R"{z) 
+ zYR"{zf -{N- l)zR'{l)R'{z) 



+ z2)i?'(l)i?"(z)) 



(27) 



where C — l/(87r^). The above expression is valid for 
both the RF and the RA models. The only difference is 
the additional inversion symmetry present in the latter: 
z goes from —1 to -t-1 in all cases, but in the RA model, 
R{~z) = R{z). 

At the same 1-loop level, the critical exponents 7y, 77 
and 7/2,772 defined in Eqs. (|20I21|) and (|25I26|) are given 



e-'Pn„n) = -(-l+6i?'(l))i?"(l) 



+ iN + 7)R"{lf +R'il)^ 



(34) 



If R{z) is analytic around z = 1, the beta functions for 
the higher derivatives evaluated at z = 1 can be derived 
as well. As noted by FisherjiS the expression for the 
pth derivative only involves derivatives of lower or equal 
order (and for p > 3 the beta function is linear in the pth 
derivative). This structure allows an iterative solution 
of the fixed-point equation, provided of course that R{z) 
has the required analytic property. 

The fixed points corresponding to the zeros of Eq. lIH^ 
are E'Jl) = (stable) and ^^1) = 1/(A^ - 2) (unsta- 
ble with an eigenvalue Ai = e). The latter fixed point 
leads to the dimensional-reduction value of the critical 
exponents, e.g., 77 = 77 = £/{N — 2), 1/ — 1/Ai = 1/e. 
The second expression, Eq. (|34|l . has then (two) nontriv- 
ial zeros only if N > 18: the fixed point with i?"(l) = 

2(jv_2)(w+7) IS unstable (A2 = 



(JV-18) X 
(Af-2) 



whereas that with i?"(l) 



(Af-8)-^(Af-2)(Af-18) 



2{N~2){N+7) 

negative second eigenvalue, A2 = — 1 



has a 



(Af-18) 



(N-2) 



This little exercise already shows that no nontrivial 
fixed-point function R^,(z), twice differentiable in z = 1, 
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can exist for A'^ < 18. Actually, one finds that there is a 
finite range of initial conditions for R'{1) for which, no 
matter what one chooses for its initial value, the RG flow 
for leads to a divergence at a finite scale /i. The 

solution to this problem has been known for some timesS^ 
the proper fixed point controlling the critical behavior 
must be nonanalytic around z = 1, with R'^{z) having a 
cusp, i.e., a term proportional to \/l — z when z ^ 1, at 
least when < 18. 

We now consider, separately and in more detail, the 
results for the RFO(iV)M and the RAO(7V)M. 



B. RF 0(7V) model 

Numerical solutions of the fixed-point equation, 
/3^, (z) = 0, have been given by FeldmanSSi for iV = 3, 4, 5 
and by us for general values of 7Vi2^ Some analytical re- 
sults can also be derived and will be discussed at the 
end of this subsection. The picture one gets from the 
numerical solutions is that the long-distance physics of 
the RF 0{N) model near d — 4 drastically depends on 
whether N is above or below two distinct critical values: 
Ndr = 18 and Nc = 2.8347.... 

The value Nun = 18 separates a region in which R'^,{z) 
at the critical, i.e., once unstable, fixed point has a cusp 
{N < Ndr) from a region {N > Ndr) where R'^,{z) has 
only a weaker nonanalyticity, a 'subcusp' in (1 — z)"^^) 
with a{N) a noninteger strictly larger than 1,22, As al- 
ready mentioned, the occurence of a cusp changes the 
values of r] and fj from their dimensional-reduction value, 
Vdr = Vdr = — 2). On the other hand, the 

weaker nonanalyticity occuring for > 18 does not alter 
the fixed-point value of i?^l) from that obtained from 
Eq. (|33|l : this leads to rj — fj — i]dr. This is illustrated in 
Figure 121 where we plot tidr/ti and tidr/ti as a function 
of N . On the other hand, at Nc = 2.834..., the pcrtur- 



riDRlri 0-6 




: 18 



N 



FIG. 2: Ratios rjon/r] (upper curve) and rjDR/fJ (lower curve) 
vs N for the RFO(N>2)M at first order in e = d - 4 
(the dimensional-reduction value for the exponents is r]DR = 
f]DR = ^/{N — 2)). The critical value of N at which both rj 
and fj diverge is A^c = 2.8347.... 



bative 'cuspy' fixed point describing the paramagnetic to 
ferromagnetic critical point when e > disappears {rj 
and fj diverge as A^ ^ A^^, see Figure Below Nc an 
attractive cuspy fixed point appears for e < that now 
describes a whole phase with QLRO. The exponents rj 
and fj characterizing this QLRO phase are plotted versus 
A^ in Figure 13 

At the critical value Nc, the beta function for R'{z) 
(unsealed by e) in exactly d — 4: has a (cuspy) fixed-point 
solution R'^{z) for any arbitrary value of the renormalized 
disorder strength R'^,{1). We have noted in Ref. [2ll that 
the situation bears some similarity with the pure 0{N) 
model near d = 2. There, the critical value Nc below 
which a QLRO phase may occur for e < is A'c = 2, 
and for Nc ~ 2 and d = 2 the beta function for the 
temperature identically vanishes, independently of the 
value of the temperature. The singular point {Nc — 2, 
d = 2) is characterized by the existence of a Kosterlitz- 
Thouless transition. One may then wonder whether the 
singular point of the RFO(Ar)M (A^^ = 2.8347..., d = 4), 
despite the absence of the Abelian property specific to the 
0{N — 2) model, also possesses a Kosterlitz-Thouless 
transition. This point will be adressed below with the 
help of the 2-loop calculation. 

We now complement this numerical study by providing 
some analytical results. We first show that for d > 4, the 
critical point is always characterized by a correlation- 
length exponent v which is equal (at one loop) to its 
dimensional-reduction value, vdr = 1/e- The eigen- 
value equation obtained by linearizing the beta function, 
Eq. (|32|l , for a small deviation from the fixed-point solu- 
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2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 

N 

FIG. 3: Exponents rj and fj (divided by e) characterizing the 
power-law decay of the pair correlations in the QLRO phase 
of the RFO{N)M for N < Nc ^ 2.8437... (first order in e = 
d-4). 
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tion, S{z) = R'{z) - R'^z), is 

^5{z) - 5{z)(r:\z)z^ + 2i?t(z)z - R'l'{z)z 
+ NR',{1) - 3i?:(l) + (4z2 + TV - 3) K'(z) - i; 
+ 5{\)({N - 3)K(z) - (TV + l)zR'l{z)- 
{z^ - 1) + 5'{z)(r:{z)z^ - 2R':'{z)z' (35) 

- NK{1)z - + 6 _ 1) R':{z)z+ 

(4z2 + iv - 3)i?:(z) + - (1 - 

By substituting S{z) = R'^{z) in the above equation, one 
can easily check that the fixed-point solution is also solu- 
tion of the eigenvalue equation with a positive eigenvalue 
Ai = e (from which, v — vdr — This result is 

independent of the analytic or nonanalytic character of 
R'^{z). For d < 4, R'^{z) is also solution of Eq. with 
Ai = e, but the eigenvalue is now negative, which allows 
the fixed point to be fully attractive. 

For iV > 18 it is possible to adapt Fisher's arguments 
concerning the hierarchy of flow equations for the succes- 
sive derivatives oi R^{z) evaluated at z = l>i^ (In his ar- 
ticle however, Fisher did not envisage nonanalytic fixed- 
point solutions.) As explained above, a fixed point with a 
well defined second derivative and an associated negative 
eigenvalue (A2 < 0) can be found for N > 18. (Note that 
cuspy fixed points are also present, but they are more 
than once unstable and correspond to putative multicrit- 
ical behavior; a detailed analysis of the fixed points and 
their stability in the N ^ 00 limit has been recently pro- 
vided by Sakamoto et aliSl) The flow equations for the 
higher derivatives are linear, namely, 

~PLd,R^^\i) = - Ap(i?'(i), i?"(i)) i?(^')(i) 

+ T,{R'{l),R"{l),--- ,R^^-'\l)), 

(36) 

provided of course that the pth derivative is well de- 
fined in z = 1. If i?'(l) and i?"(l) are chosen equal 
to their fixed-point values, i?l(l) — 1/(TV — 2) and 

^"(1) = one finds 



iV-- 2 



2p2 -{N- l)p +{N- 2)+ 
6p), 



2(A^ + 7) 



-(TV - 8 - ^/iN -2){N -18)) 



(37) 



For a given N, Ap^, monotonically increases with p, 
so that there exists an integer value p^{N) such that 
Apij(jv)* < and ^p^(N)+i* > 0- Starting with an ana- 
lytic bare action, a fixed point value is reached (provided 
R'{1) is appropriately tuned) at which the first pj(A^) 



derivatives of R^:{z) are well defined in z = 1. The RG 
flow for the (pj(TV) -I- l)th derivative on the other hand 
goes to infinity, but only in the limit ^ ^ Q. (This is 
to be contrasted with the situation for < 18 in which 
the second derivative R"{\) diverges at a finite scale /i, 
due to the nonlinear nature of the corresponding beta 
function.) As a consequence, the {p^{N) + l)th deriva- 
tive of R^{z) is not defined in z = 1, and there must be 
a nonanalyticity in R'^{z) of the form (1 — z)"^^) with 
p^{N)-l<a{N)<p^{N). 

It is easy to show that the beta function for the co- 
efficient, say a, of the (1 — z)"*^^^ term is equal to 
f3a = Ka{N)+i*a, where Ka{N)+i* is given by Eq. 
with p replaced by the noninteger a{N) -t- 1. The beta 
function is equal to zero with a nontrivial a 7^ if and 
only if Aq(jv)+i* = 0. This selects the form of the non- 
analyticity of the fixed-point solution around z = 1. As 
noticed in our previous workjSS the non-analyticity goes 
as A^/2-|-0(l) at large TV (a behavior that cannot be cap- 
tured in a 1/TV expansion, see Ref . [2^ . However, 
contrary to what stated in Ref. 22, the exponent of the 
subcusp increases continuously with N when TV > 18, as 
illustrated in Figure^ 




N 



FIG. 4: Exponent Qf(TV) characterizing the nonanalyticity (1 — 
z)"(^) in R'{z) for the RFO(TV)M in d = 4 + e (at 1 loop). 
For TV < 18 a = 1/2 whereas a(TV 18+ = 3/2. 



Finally, we close this survey of the RFO(TV)M at 1-loop 
order by pointing out that for N — 2, the FRG equa- 
tions, Eqs. (|27I29() . exactly reduce to those of a periodic 
elastic system, with a one-component displacement field, 
pinned by disorder. This is more easily seen by switching 
variables from z to an angle = cos~^(z). A cuspy fixed- 
point solution (with ry = |e|7r^/9 and fj = |e|(l -t- 7r^/9)) 
describing a QLRO phase (a 'Bragg glass') when d < 4 
is then analytically obtained (compare with Refs. I ^9l5|) . 



C. RA 0(TV) model 

Overall, the picture of the long-distance physics of the 
RA 0(TV) model that one gets from solving the 1-loop 
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FRG equations is very similar to that of the RF 0{N) 
model. This conclusion, we recall, comes with the pro- 
viso that one focuses on weak disorder (working near the 
lower critical dimension of the paramagnetic to ferromag- 
netic transition) and that the possible spin-glass ordering 
which may occur at stronger (finite) disorder is not con- 
sidered. The main difference with the behavior of the RF 
0{N) model then lies in the critical values Nun and Nc- 
Nc is found to be 9.4412...; on the other hand, Ndr = oo, 
which means that, contrary to the RF case, a cusp ap- 
pears for all values of N and dimensional reduction al- 
ways breaks down completely. 

To show that the fixed-point solution R'^,{z) has always 
a cusp, it is instructive to go back to the beta functions 
for the first two derivatives ^'(1) and ^"(1), Eqs. (|33I34|I . 
assuming that there is no cusp, (1 — z)" with < a < 1, 
in R'[z). The RA model having the additional inver- 
sion symmetry, R{—z) — R{z), it is convenient to rewrite 
R{z) = (l/2)5(z2). From Eqs. H33I34() one obtains the 
flow equations for S"(l) and S"'(l). A nontrivial fixed 
point of the beta function for S"(l) = -R'(l) is again 
i?J,(l) — 1/{N ~2) and the associated eigenvalue is posi- 
tive (Ai = e), so that the fixed point can only be reached 
if one tunes the initial value to be exactly l/(iV — 2). 
When doing so, the flow equation for S"'(l) can now be 
written as 



2(Ar-2)2 
2{N + 7)(5"(1) - Sl){S"{l) - Si) 



where 5*+ 



(jV+22)=FV(A^~2)(Af-18) 



are both strictly 



2(Af-2)(Af+7) 

negative for any finite value of N. If one starts with a 
value of >S"'(1) that is positive or even zero, which is in- 
deed a physical requirement since the bare disorder corre- 
lator is of the form A2z'^ plus possible higher-order even 
powers of z associated with even-rank anisotropics, the 
beta function, i.e., the right-hand side of Eq. H38|l . stays 
strictly positive. As a result, S"{1) diverges, and it ac- 
tually diverges at a finite scale /i. This is of course in- 
compatible with the hypothesis that the fixed point has a 
well defined second derivative S'"(l): a cusp must appear 
along the flow. Note that this reasoning is completely in- 
dependent of N. 

The limit — > oo is however somewhat special. Look- 
ing for S(z'^) — 0{1/N) and taking the N ^ oo limit in 
Eq. one finds 



2(A5"(l))oo 

((^^"(l))c 



(39) 



A fixed point with ^"(l) = 0, although having a sec- 
ond positive eigenvalue, can still be reached from an inital 
condition with S"{1) — 0: this analytic fixed point is the 
one found by a direct analysis of the N ^ oo saddle- 
point equation of an RA model with only a second-rank 



anisotropy, A2z'^. It is, however, unstable to the intro- 
duction of higher-order anisotropics (and of course never 
stable when is large but finite). 

One can actually find the cuspy fixed-point solution 
in the large A^ limit. It is convenient to change variable 
from z to (/) = cos-\z) and define f{(p) ^ {N - 2)R{z). 
T((j)) must be an even function of (f) and the symmetry 
R{—z) — R{z) translates into T{Tr — (j>) = T{(j>), so that it 
is sufficient to consider (j> in the interval [0, The beta 
function for T'((j)) in the large N limit is given by 



e-Vf,(0) = -T'(^) ^ 



cos(^^, 2. 
sin3(,/,) 



sin(0) 



(cos(20)T"(O)-f T"((/.) 



(40) 



T"(0)r"((/)) + (t"(0) - T"(0)) 



where the last term can be dropped in the large A^ limit. 
Details on the derivation of the solution and on the sta- 
bility analysis are provided in Appendix A. Here, we only 
quote the result: 



-3 sin cos 



TT+\^\ 



o 



(41) 



The 1/A^ correction can also be analytically obtained 
and is given in Appendix 1X1 One can see that Tl{(j>) in 
Eq. (|41|) satisfies the symmetry requirement around ^, 

since 3^(§) = 0, and has a cusp in \(j>\ (i.e., in \/l — z) 
as (j) goes to zero (and z goes to 1). The fixed point is 
once unstable with, as shown in the previous subsection, 
Ai = e; hence v ^ i'dr = 1/e- The critical exponents 77, 77 
and 772, ?72 are obtained from ^^1) = -f,"(0)/(A^ - 2) = 

3 I 26 
2N TP 



0(1/A3) (see Eqs. (I28l29l3()l3lll ^: 



1 



3e 
27V 
e 
2 

3e 

y 

= 2e ( 1 



(42) 
(43) 
(44) 
(45) 



Note that a Schwartz-Soffer-like inequality is satisfied, 
as it should be (see IIII by the exponents 772 and 772 
(namely, 772 < 2772), but not by the exponents rj and fj ! 

We display in Figure HI the large A^ cuspy fixed-point 
solution, Eq. H41|) . together with the numerical solution 
obtained for a wide span of N. The convergence to the 
N 00 limit is clearly visible, as is visible the presence 
for all A^'s of a nonzero slope as ^ 0, which corresponds 
to a cusp. In Figure we display the exponents 77 and 77 
as a function of A^; this again illustrates that dimensional 
reduction fails for all values of A^. 
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02 0.4 0.6 0.8 



1.4 1.6 



FIG. 5: Variation with N of the (cuspy) fixed point of the 
RAO(7V)M in d = 4 + e (at 1 loop): f."(<^) versus (p for 
< (j) < n/2. The curves correspond to A*' = 16 (top curve in 
the right part of the figure) to N — 1500 (bottom curve in the 
right part of the figure) by steps of 5% relative increase. The 
thick curve is the analytical result for N ^ oo (see Eq. 14H 'l. 
Note that there is a zero slope when <j) = 11/2 and a nonzero 
one (i.e., a cusp) when <^ ^ 0+. 




FIG. 7: Exponents r; and r] (divided by e) characterizing the 
power-law decay of the pair correlations in the QLRO phase 
of the RAO(Af)M for N < Nc ^ 9.4412... (first order in e = 
d-4). 



IV. DERIVATION OF THE FRG BETA 
FUNCTIONS TO 2 LOOPS 



4 - X 

X 

3 - X 
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N 

FIG. 6: Exponents 77 and 7? (divided by e) for the RAO(A'')M 
at first order in e = d — 4. Dimensional reduction fails for all 
values of A^: e"\ ~ 3/(2A) and e"^?? ~ 1/2 when A ^ 00. 
The critical value of N at which both 77 and 77 diverge is 
Nc = 9.4412.... 



The change from ferromagnetic ordering to QLRO oc- 
curs for Nc = 9.441..., to be compared to Nc = 2.834... 
for the RFO(Af)M. As first shown by Feldman^ QLRO 
may thus be present in the RAO(A^)M near, but below, 
d = 4 for AT = 2, 3, ...,9: see also Figure |7| The (cuspy) 
fixed-point solution associated with QLRO can be ana- 
lytically derived in the N — 2 (XY) case. Just like in the 
RFO(A^)M, the FRG equations then reduce to those of a 
disordered periodic elastic system)^*2*S the only difference 
with the RF case being a simple rescaling of the solution 
accounting for the difference in the periodicity, from 27r 
to TT. This yelds the critical exponents r/ = |e|7r^/36 and 
fj= |e|(l-K7rV36). 



In this section we describe in detail the calculation of 
the beta functions and of the critical exponents at two 
loops. 

A. Diagrammatic representation 

As explained above (see section lll Bp the calculation is 
based on an expansion of the effective action in powers 
of R. The terms of this expansion are given by all ampu- 
tated 1-particle irreducible Feynman diagrams. In order 
to determine all the counterterms, we need to compute 
the 1-replica 2-point proper vertex and the 2-replica effec- 
tive action with no derivatives (uniform fields) . The asso- 
ciated diagrams are obtained by connecting the different 
vertices of the theory with the free propagator given in 
Eq. f]l|l . 

The free propagator is represented by a line. The ver- 
tices are obtained by deriving either iSi (Eq. ITU)! or ^2 
(Eq. Illfl . In the former case they are represented by 
lines emerging from a single circle (see Fig. ^ and in the 
latter by lines emerging from two circles (corresponding 
to the two replicas of ^2) connected by a dashed line (see 
Fig. 

As can be seen in Eqs. (|8I9|) the 1- and 2-replica parts 
of the action have a factor T"^ and respectively, so 
that the various diagrams do not come with the same 
power of the temperature. Anticipating that the fixed 
point of interest to us is at zero temperature, we com- 
pute here only the diagrams of lowest order in T, i.e., 
those proportional to 1/T for the 1-replica effective ac- 
tion and to 1/T^ for the 2-replica effective action. It is 
easy to check that the diagrams of lowest order in tem- 
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FIG. 11: 1-loop diagrams for the 1-replica 2-point proper ver- 
tex. 



FIG. 8: Vertices with 3 and 4 legs obtained from Si, i.e.. 
1 rephca. 



with 



FIG. 9: Vertices with 2 and 3 legs obtained from 52, i.e.. 
2 replicas. 



with 



perature with ni 1-replica vertices and n2 2-replica ver- 
tices have rii + 2n2 — 2 propagators. Similarly, for the 
2-replica effective action, the diagrams of lowest order in 
temperature have ni + 2n2 — 1 propagators. Given these 
constraints, one can draw all the diagrams and check that 
an expansion in powers of R corresponds to an expansion 
in increasing number of loops. The 1-loop diagrams are 
displayed in Figure [TT)I and Figure ITTI 



FIG. 10: 1-loop diagrams for the 2-replica effective action. 

The 2-loop diagrams are given in Appendix ^ The 
integrals involved in the 2-loop calculations have been 
evaluated in the dimensional regularization scheme by 
using the procedure described in Appendix lUl 

The 1-replica 2-point proper vertex can be formally 
expressed as 



St:^ {q)5'n:i [q) 



H 



S2 



(46) 



1-loop -I- 2-loops 



with S = 1 — n^. The first term in the parentheses cor- 
responds to the tree diagram and is given by the inverse 
of the free propagator (see Eq. ifT^ ). 

The same can be done for the 2-replica effective action, 
which formally gives 



r2[n„ Hfa] = ^i?o(na • + 



1-loop + 2-loops. 



(47) 



We next replace the bare quantities by the renormal- 
ized ones, following Eqs. (|12|l . in the two previous expres- 
sions and reexpand in powers of R. One must then choose 



the couterterms such that the expressions are finite, i.e., 
such that all terms in 1/e and 1/e^ vanishr-- 

Once the counterterms are known, the beta functions 
can be calculated as derivatives of the renormalized quan- 
tities with respect to the scale /i at fixed bare quantities 
(see Eqs. (I16I17I18II '). In order to perform the derivative 
in Eq. (|18|l . we write the counterterm for R as: 



Zji = R + 5i[R, R] + S2[R, R, R] 



(48) 



where 5i [R, R] and S2 [R, R, R] are the 1-loop and 2- 
loop contributions, respectively; the former is a quadratic 
functional of R and is proportional to 1/e, and the latter 
is cubic in R and contains terms in 1 /e^ and 1 /e. We then 
invert Eq. H12dp with Zji given by Eq. so that the 

renormalized function R can be expressed as a functional 
of the bare function Rq as 



R = fi'^Ro—fi^'^Si[RQ, Rq] — fi^'^(S2[RQ, Ro, Ro] 
-26i[Si[Ro,RolRo])+OiR^o)- 



(49) 



The term on the last line corresponds to a repeated 1- 
loop term, in which the first argument of the functional Si 
is replaced by Si[Ro, Rq]. This gives a cubic contribution 
in Rq. The flow equation for R is then given by Eq. H18|l 
with the following beta function: 



Pr = -e{R-6i[R,R]-2S2[R, R,R] 



2Si[Si[R,R],R]). 

(50) 



In the above expression, the terms in 1 /e appearing in 
Si, S2, and in the repeated 1-loop term arc multiplied by 
e and therefore give finite contributions. The last two 
terms also have contributions in 1/e^, but they exactly 
cancel each other (as can be explicitly checked by using 
Eqs. (inSIl and below). 



B. Analytic piece of the beta functions 



As discussed in section lllll even when starting from an 
analytic initial condition, integration of the 1-loop flow 
equation generates a nonanalyticity in the renormalized 
disorder function R{z). It is therefore necessary to con- 
sider the flow equations for nonanalytic functions R{z), 
and in particular to determine the contributions of such 
nonanalyticities to the beta functions. At the 1-loop or- 
der, there are no such contributions. However, they ap- 
pear at 2 loops. For clarity's sake, we flrst describe the 
calculation of the flow equations for analytic functions, 
and consider the contributions due to nonanaliticities in 
the next section. 
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The analysis described in the previous subsection and for the counterterms: 
in Appendices B and C leads to the following expressions 



Zl, = l + ^{N- l)i?'(l) [l + ^i?'(l)(27V - 3) ) (51) 



ZJr^l + jiN - 2)i?'(l) (^1 + ^i?'(l)(iV - 2 + |) j 



(52) 



^6i[R, R]^{z^ + N~ 2) R'{z f + 2z ((1 - N)R\l) + (z^ - l) i?"(z)) i?'(z) + 4(7V - 2)i?(0)i?'(l) 
+ - 1) R!'{z) _ 1) i?"(z) _ 2i?'(l)) 



^<55[i?, i?, i?] = (e + 'i)R!\z)R!'\zf{7? - \f + 2i?"(z)2(3z(e + 4)i?"'(z) + (z^ - \)R!"\z)){z' - 1)^+ 

((9e + 32)z2 + iV(e + 2) - 2(e + %))R!\zf{z^ - \f + i?'(z)(z(15(e + 4)2^ + iV(3e + 10) - 2(6e + 25))i?"(z)2+ 

2(z2 - l)((4(e + h)7? + 2N-e- 6)R'"{z) + 2z{z^ - l)R"" {z))R" (z) + z{z^ - lf{e + 2)R"'(zf){z^ - l) + 

z((e + 4)z2 + 2N -e~ Q)R'{zf + AiN - 2)(37V + e - Q)R{z)R'{lf + 2i?'(l)2((-4A^2 _^ ^^g^ _ 9)zi?'(z) + 

(TV^^^ ^ _ 3^(^2 _ ^) _ lQ)R"(z) + (z2 - l)(2(iV + l)zR"'{z) + (z^ _ l)i?""(z))) + R'{zf{{{le + 32)^^ ^^^^ 

- 12(e + 4)z2 + 2iV2 + 5e + 2A^(z2 - l)(e + 5) + 14)i?"(z) + 2z{z'^ - l)(((e + 8)2^ + 2N -e- 6)i?"'(z) + 

z(z2 _ l)i?""(z))) - R'{l){{e + 2)R"'{zf{z^ - if + 2R"{z){z{2N + 3e + 14)i?"'(z) + 2(2^ - l)R""{z)) 

{z^ - 1)2 + ((7e + 34)2^ - 4(e + 6) + A^((e + 6)z^ + 2(e + 2)))i?"(z)2(z2 _ 1) + (_47v2 + ((e - 2)z^ - e+ 

18)iV - z^{e - 10) + e - 22)R'{zf + 2i?'(z)(z(2iV2 + (z^ - l)(e + 2)7V - e + z^{e + 18) - 20)i?"(z)+ 

(z^ _ l)(((e + I0)z2 + 2N{z'' + !)-£- 6)i?"'(z) + 2z{z'' - l)R""{z)))) 



where the superscript r stands for 'regular' i.e., analytic), 
= 87r2 and z = Tia ■ T^b + \/l — 7r^-\/l — tt^. 
Note that the expressions for Zn and correspond 
to those obtained at two loops in the nonlinear sigma 
model for the ferromagnetic-paramagnetic transition of 
the 0{N) model with no randomness if one replaces i?'(l) 
by the temperature, C by l/27r, and e by d — 2*22 More- 
over the counterterm for i?'(l) (assuming here that the 
function R{z) is analytic) reads 



-^a,(<5i+<55)U=i = ^i?'(l)(iV-2) 
1 + ^i?'(l) (iV - 2 + I 



(55) 



which again coincides with the couterterm found for the 
temperature in the pure systemiS This equivalence be- 
tween the pure model near d = 2 and the disordered one 
near d = 4 is the expression of the dimensional-reduction 
property. 



The regular part of the beta function for R{z) can then 
be calculated from Eqs. H50I53I54(I . which gives 
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= -eR(z) + ^ [(AT - 2 + z'')R'{zf - 2z((l - z'')R"{z) + {N - l)i?'(l))i?'(z) + A{N - 2)R{z)R'{l)+ 



(1 - z^)R"{z){{l - z^)R"(z) + 2i?'(l)) 



^2, 

2 



(1 - Z^) ((1 - + _ 

((1 - z2) - izR"{z) - R'{z)f + (A^ - 2)( (1 - z^fR'^zf - (l - z^) (3zi?'(z) - (z^ + 2) i?'(l)) 

- 2 (1 - i?'(z) (i?'(z) - zi?'(l)) R"{z) + (1 - i?'(l)i?'(z)2 + 4i?(z)i?'(l)2y . 



(56) 



We next consider the derivation of the critical expo- 
nents. The determination of and C,t simphfies if 
one uses the fact that Zn and Zt depend on R{z) only 
through i?'(l). Eqs. (|16I17|I become 

Ca = M^plogZ^ = logZ^ (57) 

with A being 11 or T and the derivatives being taken at 
fixed bare quantities. We then get 



C{^^{N ~ 1)CR'{1) + 0{R^) (58) 
= (iV - 2)Ci?'(l)(l + Ci?'(l)) + 0{R^). (59) 

The critical exponents ry and fj can now be evaluated 
by making use of Eqs. (|22I23|I : 

r^' ^CK{l)(l-{N -2)CK{\)) (60) 
ff =-e+{N-l)CK{l). (61) 

As we shall explicitly show in section V, the exponents 
defined under the assumption of an analytic fixed-point 
solution R<f{z) are equal to their dimensional-reduction 
value. 



C. 'Anomalous' contributions 

We have seen in section III that, even with an analytic 
initial condition for Rt=o{z), the 1-loop RG flow equation 
generates a nonanaliticity in the renormalized disorder 
function R{z). The strongest nonanaliticity is obtained 
in the RFO(A^)M for iV < 18 and in the RAO(iV)M for 
all values of N in the form 



R{z) = R{1) + R'{\){z - 1) - - (2(1 - z))* + . . . (62) 

when z approaches 1 (from below); a will be used in 
the following to quantify the strength of the singularity. 
Dimensional reduction is recovered under the assumption 
that the function i?(z) is analytic (see above), but if a 7^ 
it is no longer valid. 

Alternatively, the renormalized disorder function can 
be parametrized in terms of the angle 4> between the two 



replicas instead of the scalar product z [z — coscf}). The 
expansion in Eq. H62I) then translates into a small (j) ex- 
pansion. 



R{(j)) = R{z 



(63) 



where the nonanalyticity appears as a discontinuity in 
the third derivative of R{4>) in = 0. 

One can easily convince oneself that the nonanalytic 
term in a can explicitly appear in (3r{z). Consider the re- 
peated 1-loop term (last term in Eq. (|50|l 'l. Since 5i [R, R] 
has an explicit dependence on i?'(l), one has to compute 
dz5i\z=i which, when the nonanalyticity of R{z) is taken 
into account (see Eq. (|62|) '). takes the form 

dAU^i = f {R'{lf{N - 2) - a\N + 2)) . 

Replacing then -R'(l) by Eq.|H21 yields an explicit depen- 
dence of Pr on a. 

Actually, the 2-loop, 2-replica diagrams also give con- 
tributions in a?: look for instance at the sixth diagram 
of Figure El in the Appendix IbI Note that the two repli- 
cas of the vertices on the left and on the right of the 
diagram are actually connected via propagators. These 
vertices are therefore to be computed for identical repli- 
cas. In order to take into account the nonanalytic part 
of this diagram, the 2-replica vertices are evaluated for 
two slightly different replicas and S^/ with 



Sa' 



(64) 



Here, a is a small parameter that must be taken to zero at 
the end of the calculation and A is a unit vector orthogo- 
nal to Sa that gives the direction in which Sa' approaches 
Sa. The dependence of the diagram on A appears only 
through the scalar product A • Sf, whose absolute value 
varies between (when A and Sh are orthogonal) and 
vT^-z^ (when A, Sa and Sb are in the same plane). We 
therefore write 



A-S, 



(65) 
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with 7 varying between —1 and 1. 

There are six 2-replica diagrams giving nonanalytic 
contributions in a^: diagrams 5, 6, 11, 12, 13 and 14 of 
Figure IT^ These diagrams (and some others) also have 
contributions hnear in a, but we discard them for symme- 
try reasons. Indeed, a corresponds to the third derivative 
of R{(t>) (see Ea. (l63ll ) which changes sign under the oper- 
ation (j) —(p. On the other hand, the disorder function 
itself is unchanged in the same operation; as a result, 
linear contributions in a must vanish from all physical 
quantities. 

The situation is even more complex when considering 
the 2-loop diagrams for the 1-replica 2-point function. In 
this case a generic diagram has a singular limit a ^ 0. 
Indeed an expansion in powers of a around gives terms 
in a^a~^, a^a~^ and a^a". On top of this, there is a 
strong dependence of the result on the way the regu- 
larization is performed. Look for instance at the first 
diagram of Figure El We can decide to attribute the 
field Sa to the replica on the left of the diagram and the 
field Sa' to the replica on the right. It remains to choose 
whether the propagator (which is diagonal in replica in- 
dices) on the top of the diagram is associated with Sq or 
Sa'- The calculation shows that the two choices lead to 
different results ! 

Such ambiguities are already present in the 2-loop 
FRG treatment of disordered elastic systems and Le 
Doussal, Wiese and coworkers have given a detailed and 
well-argumented analysis of the way to handle these 
ambiguities.^^ Here, we have extended their procedure 
and used the following set of rules: 

1. The nonanalytic parts of the diagrams, i.e., those 
proportional to , come with an a priori unknown 
weight. 

2. Within a single diagram, the parts in a^a~^, a^a~^ 
and a^a'^ come with independent (a priori un- 
known) weights. 

3. For diagrams that are ambiguous in the sense that 
different regularization schemes lead to different re- 
sults, we have introduced additional weighting fac- 
tors such that all possible results can be reproduced 
by appropriately choosing these extra weighting 
factors. 

As discussed in Ref. '25, the fact that pieces of the 
2-loop diagrams come with a priori unknown weighting 
factors is due to the intrinsic ambiguity that occurs at 
T = when the function R entering into the vertices is 
nonanalytic in z = 1 (or = 0). Consider for instance 
the third derivative of R{(j>) around = (see Eg. I63|l . 
Its sign depends on whether ^ 0+ or — > 0^ and 
because of the discontinuity its value in exactly = 
is left undetermined. Vertices that contain this deriva- 
tive evaluated exactly in = have thus a contribution 
that come with an undetermined weight. Additional con- 
straints must be used to fix the values of the weighting 



factors (or at least enough relations between these fac- 
tors). This is precisely what is done by requiring that 
the physical quantities be finite. 

Under these hypotheses, the calculation proceeds in 
a straightforward way. We observe that it is possible 
to choose the nonanalytic part of the counterterms and 
to fix all the weighting factors such that the 1-replica 2- 
point proper vertex and the 2-replica effective action (see 
Eas. H46|l and H47() ') are finite. This procedure leads to a 
unique form for the counterterms (for a given parameter 
7, see Eq. (1^ ). On top of the analytic parts already 
computed (see Eas. (l51ll . (|51|l and (|54|) '). one now must 
add the 'anomalous' contributions, so that 



, 2 fiN -1){N + 2) 3N-2 



iV2 



N -2 



(66) 



(67) 



62[R,R,R] = S'2[R,R,R] 



'■{ 



(4 - iV^ - (iV - 2)e) R{z) + 
1)(N + 2) + e{3N - 

{2N + 4+(2 + j^N)e) i?"(z)}. 



^i2{N-l)iN + 2) + e{3N-2))R'iz) ^^^^ 
l-z2 



The presence of 'anomalous' contributions in the coun- 
terterms induce new terms in the beta function for R{z), 
which now reads 



/3«(z)=/3Jj(z)-^{4(iV-2)i?(z)- 

(3iV - 2)zR'{z) + (1 - z'^){-i^N + 2)i?"(z)|, 



(69) 



where, we recall, a is defined through Eq. H62|) . Note 
here that there is still an explicit dependence on 7^, 
which encodes how one takes the limit Sa' — > Sa (see 
Eq. H65|l ). There is however a preferred value for 7^. 
The simplest way to see this is to compute Pr{z) with 
a nonanalytic function R of the form given in Eq. 162|l . 
There is a term proportional to a^{l — — z. If we 

choose 7^ ^ 1, the flow equation generates a supercusp, 
i.e. a stronger nonanalyticity, R{z) ~ (1 — z)^/^, than 
the one initially considered; this supercusp would itself 
generate an even stronger nonanalyticity and the the- 
ory would not be renormalizable at 2-loop level. On the 
other hand, 7^^ = 1 ensures that the procedure is consis- 
tent. The value 7^ = 1 also appears in another context: 
the repeated 1-loop term (last term of Eq. lO) can be 
interpreted as the set of 2-loops diagrams obtained by 
replacing in the 1-loop diagrams (see Figure one of 
the 2-replica vertices by the 1-loop diagrams with two 
external legs. If we compute these 2-loop diagrams with 
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their weighting factor according to the above procedure, 
we get an expression that is consistent with the repeated 
term only if we choose 7^ = 1. In the following, we thus 
fix 7^ to this value. 

It is worth stressing that Eqs. H56I69|I (with 7^ = 1) 
exactly reduce to the 2-loop FRG equation for periodic 
disordered elastic systems when N — 2,?-5 This is more 
easily checked by switching to the angle variable (j). 

In Refs. ll2lll2l the renormalization constant for the 
temperature Zt had not been derived at two loops, so 
that an incompletely determined version of the beta func- 
tion of Eg. 1691 was given. In the notations of Ref. |21| the 
unknown parameter K is now fixed to if = 1/2 whereas 
in the notations of Ref. l28l the unknown parameter is 
now fixed to 7a = l/4i^ 

Finally, the expressions of the critical exponents have 
to be modified in order to take into account the nonana- 
lytic contributions. One obtains 

V = rf-^{N + 2) (70) 



fj = f-^m~2), (71) 
where a* is the fixed-point value for the parameter a. 

V. DISCUSSION OF THE FIXED POINTS AT 
THE 2-LOOP LEVEL AND CONCLUSION 

Going from the 1-loop to the 2-loop order does not 
significantly alter the general behavior of the RF and RA 
models in d = 4-|-e. As we shall see, it allows nonetheless 
to show that no Kosterlitz-Thouless-like transition occurs 
at the special point {Nc, d = 4) and that in the vicinity 
of this point, for N < Nc and d < A, a once unstable 
fixed point appears, which describes the transition from 
the QLRO phase to the paramagnetic one. The picture 
is now fully compatible with that found in our NP-FRG 
approach and summarized in Figure ^ 

As for the 1-loop level, it is worthwhile to start by 
considering the beta functions for the first two derivatives 
of the renormalized disorder correlator i?(z) — {e/C)R{z) 
evaluated in z = 1, assuming that the second derivative 
is well defined at that point. Writing the beta functions 
as /3 = /3i -I- /32, we only give the expressions for the 2- 
loop contributions /32, the 1-loop terms /3i being given in 
Eqs. miHH : 

e~'/3fl/(i),2 = (^-2)i?'(l)3 (72) 

e"'/^ij"(i),2 = 2(5A^ + n)R"{lf + 

6(iV + 7)R'{l)R"{lf - 6(iV - b)R'{lfR"{l) (73) 
- {N -A)R'{lf. 



One can immediately see that the fixed-point solu- 
tion -R',(l) = jv^(-'^ ~ Iv^) ^'^'^ associated (posi- 
tive) eigenvalue Ai = e(l -I- 5^32) ^^^.d to dimensional 
reduction, rj — fj — r]DR — (n-2) + W^^) (after using 

Eqs. (60,61)) and = i/Dfl, = 7 - When ^'(1) is 

chosen equal to its fixed-point value, the beta function 
for the second derivative becomes a cubic polynomial in 
R"{1) which has zeros only if the associated discriminant 
is negative, i.e., if TV > 18 - ^£,27^28 

Following the same lines as for the 1-loop order, one 
can check that the RFO(A^)M at 2 loops has a critical 
fixed point with a 'subcusp' when N > Njjfi = 18— ^e, 
but that the critical fixed point of the RAO(A^)M has 
always a cusp, implying Njjf; — 00 (see Appendix fX|l . 

We next concentrate on the region N < Nc and d < 4. 
A first result is that the beta function for R'{z) (un- 
sealed by e) in d = 4 does not vanish for arbitrary val- 
ues of the disorder strength i?l(l). When scaling out 
the disorder strength to define r'{z) = R'{z)/R'^{l), the 
beta function for e = can be expressed as (3r'{z) = 
/3i[r',r'] + R'^{l)l32[r',r',r']. f3r'{z) then vanishes inde- 
pendently of the value of -R'* (1) if and only if (3i [r', r'] = 
and (32[r' ^r' ,r'] — have the same solution r'^{z). It 
is straightforward to check that inserting the solution 
of the 1-loop equation /3i = into the 2-loop equation 
does not make the latter vanish identically. The conse- 
quence is that no Kosterlitz-Thouless transition can exist 
in {Nc,d = 4) since, even perturbatively, no line of fixed 
points can be found. Actually, the only fixed point at 
2-loop level for N ~ Nc and d = 4 is the trivial one with 

Kil) = 0. 

For N ^ Nc and d < 4, consideration of the 2-loop 
order brings in a new phenomenon. An additional, once 
unstable, fixed point appears, which describes the tran- 
sition from the QLRO phase to the paramagnetic one. 
This new fixed point is found perturbatively in e (which 
is now negative) and N — Nc. More precisely, as shown 
by Le Doussal and Wiese^^ it can be obtained within 
a double expansion in ^/\e\ and N — Nc. For any given 
value N ^ Nc, the critical fixed point and the QLRO 
fixed point get closer as |e| increases and they merge for a 
value eic{N). This latter corresponds to the lower critical 
dimension d/c = 4 + eic of QLRO. With the full 2-loop 
results given above, Eqs. H56I69() . one finds 

dic(RF) = 4 - 0.14(iV - Ncf + 0{{N - Ncf) (74) 
dic(RA) = 4 - 0.002(iV - Nc f + 0{{N - Ncf) (75) 

where Nc = 2.8347... for the RF model and 9.4412... 
for the RA model. For what it is worth, directly plug- 
ging iV = 2 in the above expressions gives the fol- 
lowing estimates for the lower critical dimension of the 
(QLRO) Bragg glass phase in the XY model: dic(RF) ~ 
3.9,dic(RA) ~ 3.9. (Our NP-FRG theory predicts 
dic(RF) ~ 3.8,^^ so that in all cases no Bragg glass is 
found for iV = 2 in d = 3.) 

Note that the fact that no Kosterlitz-Thouless tran- 
sition takes place in N = Nc and d = 4 is con- 
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nected to the absence of a linear term in N — Nc in 
Eqs. H74I75|I . (In the pure 0{N) model on the contrary, 
dia^2- 6(2 - N) + 0{{N - 2)2), where 6 = 1/4 is pro- 
portional to the temperature of the Kosterlitz-Thouless 
transition.'-^) The slope of the curve Nic{d) giving the 
locus of the lower critical dimension for QLRO is infinite 
as d ^ 4~ : see Figure ^1 

The 2-loop predictions around N — Nc and N — 
Nic{d) near d = 4 are in agreement with those of our 
NP-FRG treatment. For a detailed comparison, we plot 
the two sets of results in Figures I|12I13|) . The nonper- 
turbative, but approximate FRG computation does not 
exactly reproduces the 2-loop calculation near d = 4, but 
the differences are not very significant and do not alter 
the general picture. 
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RA 0{N) models near d = 4 is perturbatively renormal- 
izable at two loops, thereby proving that the 1-loop result 
is not fortuitous. The results we have obtained within 
the 2-loop order fit into the general scenario predicted 
by our NP-FRG approach. ^^-^^ Considering the techni- 
cal difficulties associated with the FRG loop expansion 
in d = 4 -I- e, it is highly unlikely that perturbative FRG 
will ever provide accurate extrapolations to the physical 
cases d = 2, 3 (and to = 1) for the RF and RA mod- 
els. The NP-FRG on the other hand offers a direct way 
to study these situations. 



APPENDIX A: FIXED POINTS AND THEIR 
STABILITY FOR THE RAC)(iV)M IN THE LARGE 
N LIMIT 

We first rewrite the 1-loop beta function for the deriva- 
tive of f (0) = {N-2)R{z), Eq. gUJ), by introducing the 
function U{(p) = -^{(p)/ sin cj), 

sin(0)/3t;(0) ={U{0) - 1) sin(0)[/(0) + 

(f/(0)cos(</.) -[/(</.))[/'(</)), 

and we look for fixed-point solutions. When t/,(0) = 1, 
the only solutions are f/*((/)) = cos^ and /7*(0) = 1. If 
L/*(0) ^ 1, the equation (3u{4') = can be solved by 
inverting the relation between U and 4> s-nd considering 
(j) as a, function of U . The fixed-point equation now reads 



FIG. 12: Comparison between the results of the 2-loop per- 
turbative FRG (dashed line) and of the NP-FRG (continu- 
ous line) near d = 4 for the QLRO lower critical dimension 
Nic{d) in the case of the RFO(iV)M. The black circle denotes 
the physical case of the XY model in d = 3, a case which is 
clearly below its lower critical dimension. 




FIG. 13: Comparison between the results of the 2-loop per- 
turbative FRG (dashed line) and of the NP-FRG (continuous 
line) near d = 4 for Noaid) in the case of the RFO(iV)M. 

To summarize: We have shown in this paper that the 
theory describing the long-distance physics of the RF and 



5c/ (cos ^(L/)) 



C/o cos(0(t/)) 



1 



U 



C/n - 1 



(A2) 



where C/o is such that 0(f/o) — 0- The solutions of 

Up 

Eq. HA2p are easily found as cos((/)(C/)) = KU"o-^ + JJ. 
K \s & constant that is determined through the condi- 
tion 4){Uq) = 0. The result can be reexpressed by stating 
that the fixed-point functions U<t{(j)) are solutions of the 
following transcendental equation: 



c/*(o) 



cos((/)) (A3) 



where C/*(0) is different from 1 but still unknown. 

Note that the property T{Tr—(j>) = T{(t)) imposes U{tt— 
(j)) = — [/{(/)). Since we expect the functions to be analytic 
around | {i.e., z = 0), this property implies that C/(f ) 
and all even derivatives of U evaluated in ^ are equal 
to zero. This requirement can only be fulfilled by the 
solutions of Eq. (|A3I) if ^^^"J^^i is a nonzero integer. On 
the other hand, it cannot be satisfied by the solution 
U„{(j)) — 1, which should therefore be discarded. 

To determine the acceptable values of J7*(0), we have 
to turn to the stability analysis. Introducing a small 
deviation around the fixed point, S{(f>) = U{<j)) — U^:{4>), 
and linearizing the associated FRG equation leads to the 
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following eigenvalue equation: 

(e-^A) sin0 S{(j)) = (sin0 [/, (,^) + cos </> U!,{(f>))S{Q) 
+ ((C/.(O)-l)sin0-C/:(^))<5(^) 
+ (C/*(0) cos <?!>-[/,(</.) 

(A4) 

We first check that the solution U^:{4>) — coscj), corre- 
sponding to [/*(0) = 1, is fully unstable with e^^A — 1, 
and we next consider the solutions given by Eq. ^M^. As 
before we consider (p as a function oi U = C/*((/i), which 
gives after some manipulations, 



-{U, - 1) I 1 - I j US'iU) 



From this expression, one derives — T"(0) = |G'(^) = 
I + ^ + 0{-^ ), which leads to the expressions given in 
section Fill CI 

Finally, the fixed point can also be found at the 2-loop 
level in the large N limit by using the above variable 
X. One obtains for instance that the correlation length 
exponent is equal to 



1 
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37V 



29707 

81iV2 



(A9) 



so that it is now different from the dimensional-reduction 
value. 



APPENDIX B: 2-LOOPS DIAGRAMS 



1 - (e-^A - C/o + 1) 1 




S{U) 



(A5) 



We give here all the 2-loop diagrams built with the 
rules defined in section IV A. 



where Uq = U^,{(f> — Q). The solution of the above equa- 
tions is easily obtained as: 



APPENDIX C: 2-LOOP INTEGRALS IN 
DIMENSIONAL REGULARIZATION 



"0- 



S{Uo) f U\ 

F^A [iTn 



1 - 




(A6) 



where, we recall, Uq/{Uq — 1) is a nonzero integer. The 
condition that 5{U) be an odd fuction of [/, analytic 
around t/ = 0, imposes stringent constraints on e~^A. 
One finds that the only singly unstable fixed point corre- 
sponds to Uo/{Uo — l) = 3, i.e., L/*(0) = 3/2. The associ- 
ated eigenvalues are equal to e~^A = 1, 0, —1, —2, —3, .... 
For C/*(0) = 3/2, Eq. HA3p can be exphcitly solved, which 
leads to [/*((/)) = 3cos(^) and to Eq. (gTJ. 

To derive the fixed-point solution at the following or- 
ders in 1/iV, great simplification is obtained by first notic- 
ing that the above solution (in the limit N oo) can be 
rewritten as 



3 ,7r-20. 



2cos(^-^^)-l 



(A7) 



which implies that is a function of cos( ^ g '*' ). For 

finite N we now look for a fixed-point solution of the 
form fi{(t)) = _|sin(^^)(2X - l)G(X) with X = 



cos( ^~ ''' ). G(X) can be expanded in powers of 1/iV 
(or for convenience, 1/(A^ — 2)), its leading term being 
simply equal to 1 (see Eq. ljA7p ). With this transfor- 
mation the fixed-point equation can be solved in powers 
of l/(iV — 2), each term being a polynomial in X. One 
obtains for the first terms 

G{X) = 1 + (95 - AAX - IQX^) 



<d{N -2) 



O 



1 



(A8) 



In this section, we discuss the procedure used to eval- 
uate the integrals appearing in the two-loop calculation 
in d = 4 + e. 



1. Integrals for the 2-replica diagrams 

We start with the integrals for the 2-replica effective 
action. The most general integral reads 



J. 



i,j,k,l,'m,n(^^ ^) ~ 

(qi 



(5(qi +q2 +q3) 



q2)''(q2 • q3)'*(q3 • qi)* 



ill 



.ay{ql + anql+by{ql + 



b)^{ql + by 



(CI) 



with the parameters {i, j, k,l,m,n,r, s, t} being nonneg- 
ative integers, {a, 6} being positive real numbers and 
r + s + t < 2. The integrals have obvious symmetry prop- 
erties since the integration variables can be exchanged. 
For instance, the integral is unchanged when {z, I, s} and 
{j,m,t} are exchanged, when {j, m,r} and {k,n,t} are 
exchanged, or when fc, a} and {l,m,n,b} are ex- 
changed. 

It is possible to reduce the range of values of {r, s, t} 
that one must consider. This can be done by rewriting in 
the integrand qi • q2 = {qi — qf — '?2)/2- Qi can then be 
replaced by {qf + a) — a, which, if i > 0, can be combined 



with {qf + a) Ho give {qf 



-i+l 



— a{qf+a) Similar 



transformations can be done for and q^. Then, under 
the assumption that i, j and k are nonzero, one gets the 
identity 



jrst 



1 



V i,j.k — l,l,rii,n 



r — l,s,t 



i^j — l,k.l,m.n 



'l,S,t 



(C2) 



i,j,k,l,7n,nJ 



17 



T ! 1 T yi 1 t- 



l k 



f /I 1 r 





i I i I i I i 



'3 



FIG. 14: 2-loop diagrams for the 2-replica eflective action. 



FIG. 15: 2-loop diagrams of type A for the 1-replica 2-point proper vertex. 

J!5;f JZST )C3 X3 ^3 



FIG. 16: 2-loop diagrams of type B for the 1-replica 2-point proper vertex. 



FIG. 17: 2-loop diagrams of type C for the 1-replica 2-point proper vertex. 
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FIG. 18: 2-loop diagrams of type D for the 1-replica 2-point proper vertex. 

Jr3 dZS 13 JZ^^ J3- JZ!r JZS JZi C3 



izT IZT iz^F iz^P [3= [3= 

FIG. 19: 2-loop diagrams of type E for the 1-replica 2-point proper vertex. 
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FIG. 20: 2-loop diagrams of type F for the 1-replica 2-poiiit proper vertex. 



Similar relations can be obtained under the (weaker) con- 
dition that i + I, j + m and k + n are nonzero. 

There remains to treat the case where one of the three 
previous combinations (say for instance the last one) van- 
ishes (which implies k ^ n ~ 0). It is sufficient to treat 
the three cases where {r,s,t} = {1,0,0} up to a per- 
mutation, because the other possibilities never appear in 
the calculation. Consider for instance Jj^^*!, , ^ q. Since 
53 appears in the integrand only through the S function, 
one performs the integral on qs trivially. The integrand 
can then be factorized in a piece qi((Ji + a)^'^{qf + 
depending only on qi and another (of a similar form) 
depending only on q2. Each piece is a vector so that the 
integral vanishes for symmetry reasons. Consider now 
the integral J^fj^o i mO- integrating over qs, the nu- 
merator of the integrand becomes — (^2 + Qi • 12 )• The 
last term gives zero after integration for symmetry rea- 
sons, just as before. The first term can be rewritten as 
(q| -|- a) — a so that, under the condition that j > 

7OIO _ „ tOOO _ tOOO 



(C3) 



Similar equations can be obtained if j = and m > 0, or 
if r = s = 0, t = 1. 

We can further simplify the integrals by using the re- 
lation 



1 



1 



(g2 +a)(q^ + h) a-b\q^ + h + 



1 



1 



(C4) 



which enables one to reduce the integrals to a form where, 
in the three couples {«,/}, {j, {fc,n}, at least one el- 
ement is zero. 

The previous procedure reduces the problem to the 
case r = s = t = Q. In this case, the integral is finite 
(and therefore of no interest for us since we work in the 
minimal subtraction scheme) whenever the four condi- 
tions i + j + k + l + m + n > 4, i + j + I + m, > 2, 
j + k + m + n>2, k + i + n + l>2 are satisfied simul- 
taneously. The divergent integrals can be of two types: 

1. if the smallest value between i + I, j + m, fc -|- n is 

zero, then the integral factorizes, and identifies to 
a product of 1-loop integrals of the form: 



Ii{a) 



L 



1 



1 T{i-d/2) 



,d/2- 



(C5) 



,q(g2+a)' {47r)d/^ T{i) 

2. if the smallest value between i + I, j + m, k + n 
is 1, then the second smallest one must also be 1 



for the integral to be divergent, and the last one is 
free. All these integrals can be calculated by using 
derivatives of the relation: 



i 



Sjqi + q2 + qs) 

qiq2q3 



-C 



d-3 




(C6) 



with = Stt^. 



2. Integrals for the 1-replica diagrams 

The situation for the 1-replica diagrams is in a sense 
simpler because only one mass can appear. However, the 
vertices and propagators can now come with the external 
momentum p. The first step in the evaluation of the 
integrals consists in expanding the integrand in powers of 
p, keeping only the constant, linear and quadratic terms 
which are of interest in the renormalization procedure. 
The linear part of the integral vanishes for symmetry 
reasons. The constant part can be evaluated by using 
the procedure described in the previous subsection. 

The quadratic part can come in the two following 
forms: p"^ or (q^ • p)(qj • p). The first case is simple 
and can be evalated by using the procedure described 
above. For the second form, we used the relation 



qiq2q3 



(^(qi + q2 + q3)(qi • p)(qi • p)/ = 

p2 f 

— / (5(qi +q2 + q3)(qi • q?) / 

^ •'qiq2q3 



(C7) 



where / is a function of qi, q2 and qs and i,j are equal 

to 1, 2 or 3. The case i 7^ j can then be treated as in the 
previous section. In the case i = j, the integrand always 
comes with {qf + c)" in the denominator, so one can use 
again the relation: 



1 



(C8) 



and compute the remaining part by using the procedure 
described previously. 
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